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The  gas-phase  chemistry  of  the  hypergolic  system  CH3NHNH2  -  monomethylhydrazine  (MMH),  with  oxi¬ 
dizers  NO2/N2O4  at  room  temperature  and  1  atm  N2  was  investigated  experimentally  using  a  gold-coated 
chamber  reactor,  coupled  with  a  Fourier  transform  infrared  (FTIR)  spectrometer.  The  IR-active  species 
identified  in  the  early  reactions  include  HONO,  monomethylhydrazinium  nitrite  (MMHHONO),  methyl 
diazene  (CH3N=NH),  methyl  nitrate  (CH30N02),  methyl  nitrite  (CH3ONO),  nitromethane  (CH3N02), 
methyl  azide  (CH3N3),  H20,  N20  and  NO.  In  order  to  elucidate  the  mechanisms  by  which  these  observed 
products  are  formed,  we  carried  out  quantum  mechanics  calculations  [CCSD(T)/M06-2X]  for  the  possible 
reaction  pathways.  Based  on  these  studies,  we  propose  that  the  oxidation  of  MMH  in  an  atmosphere  of 
N02  occurs  via  two  mechanisms:  (1)  sequential  H-abstraction  and  HONO  formation,  and  (2)  reaction  of 
MMH  with  asymmetric  0N0N02,  leading  to  formation  of  methyl  nitrate.  These  mechanisms  successfully 
explain  all  intermediates  observed  experimentally.  We  conclude  that  the  formation  of  asymmetric 
ONONO2  is  assisted  by  an  aerosol  formed  by  HONO  and  MMH  that  provides  a  large  surface  area  for 
ONONO2  to  condense,  leading  to  the  generation  of  methyl  nitrate.  Thus  we  propose  that  the  overall 
pre-ignition  process  involves  both  gas-phase  and  aerosol-phase  reactions. 

©  2013  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 


1.  Introduction 

Hypergolic  bipropellants  are  fuel-oxidizer  combinations  that 
ignite  spontaneously  upon  mixing  at  ambient  temperatures.  They 
facilitate  the  design  of  rocket  thrusters  by  simplifying  the  ignition 
system,  and  are  widely  used  in  propulsion  systems  in  which  vari¬ 
able  and/or  intermittent  thrust  capabilities  are  needed.  Among  the 
most  commonly  deployed  bipropellant  combinations  is  mono- 
methylhydrazine/nitrogen  tetroxide,  which  is  also  referred  to  as 
MMH/NTO  or  CH3NHNH2/N2O4  [1]  For  applications  in  which  the 
freezing  point  of  NTO  is  too  high,  an  alternative  oxidizer  is  red 
fuming  nitric  acid  (RFNA),  which  is  composed  of  nitric  acid 
(HN03,  ~85  wt%)  and  N02  (8-15  wt%). 

Recently  the  impinging  stream  vortex  engine  (ISVE)  has  at¬ 
tracted  significant  attention  due  to  its  compact  size  and  potential 
for  efficient  combustion,  making  it  important  to  develop  computa¬ 
tional  fluid  dynamics  (CFD)  models  to  gain  insight  into  the  influ¬ 
ence  of  design  parameters  on  engine  performance  [2-5].  An 
important  part  of  this  effort  is  to  develop  a  chemical  kinetics  mech¬ 
anism  for  MMH/NTO  or  MMH/RFNA  combinations. 
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To  provide  a  starting  point  for  such  activities  a  detailed,  finite- 
rate,  chemical  kinetics  mechanism  of  MMH/RFNA  was  developed 
by  the  U.S.  Army  Research  Laboratory  (ARL)  [6-9]  for  modeling 
the  gas-phase  combustion  processes.  The  most  recent  version  of 
this  mechanism  involves  513  reactions  and  81  species  [6].  Sources 
for  the  ARL  mechanism  include  the  following: 

1.  a  set  of  reactions  for  H/C/N/O  compounds  developed  by  Ander¬ 
son  and  co-workers  for  modeling  the  dark  zones  observed  in 
solid-propellant  combustion  (43  species,  204  reactions)  [10], 

2.  approximately  160  small-hydrocarbon-molecule  reactions  that 
were  extracted  from  the  GR1  3.0  database  [11], 

3.  approximately  80  reactions  involving  HN03,  N03,  N204,  and 
hydrocarbon/NOx  moieties  that  were  identified  via  a  literature 
search  performed  specifically  for  the  mechanism  development 
effort,  and 

4.  approximately  50  reactions  recommended  by  Catoire  and  co¬ 
workers  for  modeling  the  ignition  and  combustion  of  MMH/ 
02  [12]  and  MMH/NTO  [13]  systems. 

The  validity  and  completeness  of  the  ARL  mechanism  was 
tested  by  running  CHEMK1N  [14]  simulations  for  MMH/NTO  sys¬ 
tems,  and  a  reduced  version  of  the  mechanism  was  used  in  CFD 
simulations  for  the  ISVE  engine  [2,5,8]. 
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One  major  concern  with  the  ARL  MMH/RFNA  mechanism  is  the 
lack  of  relevant  experimental  studies  for  its  validation  [6].  As  part 
of  an  effort  to  provide  experimental  support  for  this  mechanism, 
we  investigated  the  pre-ignition  reactions  between  MMH  and 
HN03  (the  major  constituent  of  RFNA)  in  an  earlier  work  [15]. 
These  experimental  results  suggested  that  the  current  MMH/RFNA 
mechanism  omits  some  important  early  reactions  between  MMH 
and  HN03  and  corresponding  species.  Since  N02/N204  is  another 
important  constituent  in  RFNA,  its  early  gas-phase  reactions  with 
MMH  are  examined  in  this  work. 

The  current  ARL  mechanism  for  MMH/NTO,  a  subset  of  the 
MMH/RFNA  mechanism,  contains  reactions  categorized  in  two 
domains: 

1.  single-bond  fission  events  to  strip  fragments  from  MMH  and 
generate  free  radicals,  and 

2.  radical-radical  reactions  to  form  either  closed-shell  or  open- 
shell  species. 

Given  the  low  temperature  (<100  °C)  in  the  pre-ignition  envi¬ 
ronment,  direct  bond  fission  from  MMH  to  produce  H,  CH3  or 
NH2  is  unlikely  so  that  N02  is  the  major  free  radical  available  ini¬ 
tially.  Based  on  this  assumption,  the  ARL  mechanism  considers 
the  two  types  of  initial  reactions: 

1.  H-abstraction  from  MMH  and  sequential  HONO  formations,  and 

2.  recombination  between  N02  and  MMH  free  radicals  generated 
by  H-abstraction. 

However,  these  reactions  do  not  fully  explain  the  formation  of  a 
condensate  that  has  been  observed  in  several  previous  studies 
involving  examinations  of  a  residue  from  gas-phase  reactions  in  a 
stoichiometric  mixture  of  MMH  and  N02/N204  [16].  The  IR  proper¬ 
ties  of  this  residue  are  quite  similar  to  the  IR  properties  of  a  residue 
obtained  from  reactions  between  liquid-phase  MMH  and  gaseous 
N02/N204  [17]  The  IR  properties  of  the  residue  from  these  two  stud¬ 
ies  suggest  that  monomethylhydrazinium  nitrate  (MMH  HN03)  is 
formed  in  addition  to  other  species  [16].  The  formation  of 
MMH  HN03  was  also  detected  by  Saad  et  al.  [18],  who  examined  li¬ 
quid-phase  reactions  between  MMH  and  N204  in  a  system  diluted 
by  CC14  at  -20  °C.  In  a  recent  work  by  Catoire  et  al.  [13],  it  is  sug¬ 
gested  that  the  MMH  HN03  detected  as  a  major  product  in  the  res¬ 
idue  by  Semans  et  al.  is  not  formed  from  reactions  in  the  gas  phase, 
since  its  elemental  analysis  matches  rather  poorly  with  that  of  the 
residue  as  determined  by  Breisacher  et  al.  [19].  Catoire  et  al.  sug¬ 
gests  that  nonionic  compounds  are  formed  and  accumulate  in  a 
condensate.  However,  it  appears  that  no  experiments  were  carried 
out  to  confirm  the  formation  of  these  nonionic  compounds. 

Based  on  the  above  discussion,  there  is  a  clear  need  to  reexam¬ 
ine  gas-phase  reactions  between  MMH  and  N02/N204  at  low  tem¬ 
peratures  in  order  to  identify  the  relevant  preignition  products  and 
reaction  pathways.  There  are  two  objectives  with  the  present  work. 
First,  we  would  like  to  experimentally  identify  in  situ  the  species 
formed  early  in  the  preignition  event  from  gas-phase  reactions  be¬ 
tween  MMH  and  N02/N204.  Second,  we  would  like  to  use  quantum 
mechanics  (QMs)  tools  to  help  elucidate  the  reaction  pathways, 
since  experimentally  it  is  rather  difficult  to  identify  and  quantify 
radicals,  as  well  as  to  identify  transition-state  structures. 


2.  Methods 

2.1.  Chemicals 

Monomethylhydrazine  (MMH),  purchased  from  Sigma- 
Aldrich  and  used  without  further  purification,  is  a  liquid  with  a 


concentration  higher  than  98%.  Nitrogen  dioxide  (N02)  is  diluted 
in  nitrogen  (N2)  and  stored  in  a  pressurized  cylinder.  The  mole  con¬ 
centration  of  N02  is  1%.  In  the  gas  phase,  nitrogen  dioxide  and  dini¬ 
trogen  tetroxide  always  coexist  due  to  the  reversible  dimerization 
reaction:  2N02  <->  N204  [20].  Since  MMH  and  N02/N204  are  highly 
toxic  and  corrosive,  only  small  amounts  of  reactants  were  used  in 
each  test  (i.e.  ~0.2  pi  MMH)  and  necessary  protections  (i.e.,  gloves, 
mask,  etc.)  were  adopted. 

2.2.  Experimental  setup  and  procedures 

An  overview  of  the  experimental  setup  is  shown  in  Fig.  1.  The 
chamber  reactor  is  made  of  stainless  steel  and  has  a  12  mm  inner 
diameter  and  a  path  length  of  65  mm.  A  thin  layer  of  gold  is  coated 
on  the  chamber  surface  to  minimize  the  potential  of  surface  cata¬ 
lytic  reactions.  On  each  side  of  the  cylindrical  chamber  is  a  ZnSe 
window  for  transmission  of  the  modulated  FT-IR  beam  from  a  Bru- 
ker  IFS  66/S  FT-IR  spectrometer.  N02/N204/N2  mixture  is  injected 
into  the  chamber  reactor  through  a  Teflon  tube,  while  MMH 
(~0.2  pi)  is  delivered  by  a  1  pi  syringe.  MMH  is  volatile  (vapor 
pressure  is  5  kPa  at  20  °C)  and  evaporates  rapidly  upon  injection. 
MMH  vapor  then  mixes  and  reacts  with  N02/N204  in  the  chamber. 
The  temporal  evolution  of  IR-active  species  in  the  chamber,  includ¬ 
ing  initial  reactants  and  products,  is  recorded  by  a  time-resolved 
rapid  scan  on  the  FT-IR  spectrometer.  In  each  test  a  total  of  390 
spectra  (the  maximum  memory  limit),  with  a  spectral  resolution 
of  2.6  cm~'  and  temporal  resolution  of  50  ms,  are  acquired.  It 
should  be  noted  that  IR-inactive  species  such  as  N2,  02  and  H2, 
transparent  to  the  IR  beam,  cannot  be  detected.  Since  the  reactants 
are  highly  diluted  in  N2  (~99%),  the  temperature  and  pressure  vari¬ 
ations  in  the  chamber  due  to  chemical  reactions  can  be  neglected. 

Two  types  of  tests,  labeled  as  either  oxidizer-rich  or  fuel-rich  in 
the  discussion  below,  were  conducted. 

•  In  the  oxidizer-rich  test,  the  reactor  is  purged  with  N02/N204/ 
N2  mixture  and  thus  initially  filled  with  the  oxidizer.  0.2  pL 
MMH  (3.8  x  10~6mol)  is  quickly  injected  into  the  chamber 
and  reacts  with  N02/N204. 

•  In  the  fuel-rich  test,  the  chamber  is  purged  with  N2  and  1  pL 
MMH  (1.9  x  10~5  mol)  is  first  injected  into  the  chamber,  which 
is  then  initially  filled  with  MMH  vapor  and  N2.  The  mixture  of 
oxidizers  (N02/N204/N2)  is  then  pressurized  into  the  chamber 
and  reacts  with  MMH  vapor  upon  mixing.  The  findings  from 
these  two  types  of  tests  are  discussed  separately  in  the  results 
section. 

Because  residues  may  remain  on  the  chamber  wall  after  each 
test,  the  chamber  wall  is  cleaned  and  conditioned  through  two 
steps  before  each  new  test.  In  the  first  step  the  chamber  reactor 
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is  heated  to  about  350  °C  by  two  embedded  cartridge  heaters.  Most 
residues  will  either  be  gasified  or  decompose  at  this  temperature, 
and  therefore  removed  by  the  purge  gas.  In  the  second  step  the 
chamber  wall  is  conditioned  with  N02/N204  or  MMH  vapor, 
depending  on  which  one  is  introduced  first.  The  conditioning  is  re¬ 
garded  complete  when  there  is  no  noticeable  decrease  of  N02/N204 
(or  MMH)  in  the  closed  chamber  after  a  certain  period  of  time. 

2.3.  Computational  details 

The  geometry  optimization  and  Hessian  calculation  were  car¬ 
ried  out  at  the  level  of  M06-2X/6-311++G**  [21].  The  Hessian  was 
used  to  provide  the  vibrational  frequencies  for  zero-point  energy 
(ZPE)  and  thermocorrections  to  enthalpy  and  entropy.  In  addition, 
at  these  optimized  geometries  we  calculated  the  energy  at  the 
UCCSD(T)/6-31G”  level  of  QM.  In  the  reaction  of  HONO  formation, 
the  ONO-H  distance  is  the  key  reaction  coordinate  and  sensitive  to 
different  functionals.  Comparing  with  geometries  at  transition 
state  (TS)  from  CCSD/6-31+G**  reported  by  McQuaid  and  Ishikawa 
[7],  the  greatest  difference  in  O— H  distance  is  at  the  TS  of  reaction 
to  from  CH3NNH2+HONO  (2.170  v.s  1.908  A),  and  for  the  remain¬ 
ing,  the  difference  in  0— H  is  less  than  0.1  A.  These  geometries  dif¬ 
ferences,  however,  do  not  cause  much  difference  in  barrier  heights 
of  HONO  formation,  as  shown  in  Section  4.1.  All  TS  were  shown  to 
have  exactly  one  negative  eigenvalue  by  following  the  minimum 
energy  path  (MEP)  scan  to  connect  reactant  and  product.  Free  ener¬ 
gies  are  reported  at  298.15  K  and  1  atm. 

For  reactions  inside  the  aerosol  or  on  the  aerosol  surface,  we 
added  the  electrostatic  interaction  between  reactants  and  sur¬ 
rounding  ions  with  Poisson-Boltzmann  solvation  model  (implicit 
solvent)  implemented  in  Jaguar  [22],  using  a  dielectric  constant 
of  80.37  and  a  spherical  cavity  of  radius  1.40A  for  water.  We  con¬ 
sider  that  the  solvation  effects  calculated  for  water  represent  the 
high  dielectric  properties  expected  for  these  systems,  with  the  re¬ 
sults  depending  little  on  the  exact  values  as  long  as  the  dielectric 
constant  is  greater  than  20  and  radius  smaller  than  2.8A.  All  geom¬ 
etry  optimizations,  solvation  and  Hessian  calculations  were  carried 
out  with  Jaguar  7.6  [23].  The  UCCSD(T)  calculation  was  done  with 
NWChern  [24,25], 

3.  Experimental  results 

3.1.  Summary  of  spectroscopic  data 

Figures  2a,  b  and  c  show  selected  FT-1R  spectra  obtained  from  a 
typical  test  with  the  chamber  initially  filled  with  the  oxidizer  mix¬ 
ture  (N02/N204/N2);  that  is,  the  oxidizer-rich  test.  As  shown  in 
Fig.  2a,  a  chemical  equilibrium  is  reached  between  N02  and  N204. 
The  equilibrium  constant  I<p  =  P2(N02)/P(N204)  is  about  0.1  atm  at 
room  temperature  (20  °C)  [26].  Therefore  the  initial  mole  composi¬ 
tion  of  N02/N204/N2  is  about  0.855%/0.073%/99.072%  at  room  tem¬ 
perature  and  1  atmosphere  in  the  chamber.  After  injecting  0.2  pi  of 
MMH,  the  concentration  of  N02/N204  decreases  and  the  products 
increase  with  time.  Figures  2b  and  c  show  spectra  acquired  at 
t  =  7.5  s  and  15  s,  respectively  (t  =  0  at  MMH  injection).  The  injec¬ 
tion  of  the  MMH  takes  a  fraction  of  a  second. 

The  original  IR  spectrum  shown  in  Fig.  2c  is  a  spectral  summa¬ 
tion  of  all  IR-active  species,  including  reactants  and  products.  Most 
species  can  be  readily  distinguished  by  their  characteristic  vibra¬ 
tional  frequencies.  The  identified  species  are  labeled  in  Fig.  2c 
and  their  vibrational  frequencies  specified  in  Table  1.  To  identify 
other  species  in  the  product  spectrum  further  spectral  analysis 
was  conducted  by  subtracting  spectral  bands  of  H20,  N02,  N204, 
N20,  and  CH30N02  with  the  remaining  spectrum  shown  in 
Fig.  2d.  Spectral  subtraction  is  an  important  method  to  simplify 


the  interpretation  and  analysis  of  the  product  spectrum  [27].  For 
comparison,  a  spectrum  of  monomethylhydrazinium  nitrate 
(MMH  HN03),  which  was  acquired  in  an  earlier  study  on  MMH/ 
HN03  reactions,  is  shown  in  Fig.  2e.  The  broad  spectral  bands  at 
3322,  3080,  2796  and  2495  cm-1  in  Fig.  2d  match  well  with  those 
in  Fig.  2e,  and  these  spectral  bands  were  assigned  to  monometh- 
lydrazinium  cation  which  may  exist  in  three  forms  [15]: 
[CH3NHNH3]+,  [CH3NH2NH2]+  and  [CH3NH2NH3]2+. 

The  spectral  bands  in  the  region  1000-1500  cm~'  do  not  match 
in  Fig.  2d  and  e.  In  Fig.  2e,  one  band  is  centered  at  1380  cirr1,  and  it 
is  typical  for  a  nitrate  salt  [28,29].  In  Fig.  2d,  this  band  is  centered 
at  1215  cm~\  where  many  inorganic  nitrite  salts  have  a  broad 
band  [28-30].  Therefore,  the  broad  spectral  bands  at  2400- 
3400  cm-1  and  1000-1500  cur1  can  be  assigned  to  mono¬ 
methylhydrazinium  nitrite  (MMH  HONO).  In  addition,  several 
spectral  bands  still  exist  in  the  region  near  1600  cur1  after  sub¬ 
tracting  bands  attributed  to  N02.  These  remaining  bands  are  pri¬ 
marily  from  nitromethane  (CH3N02)  and  methyl  nitrite  (CH3ONO). 

A  summary  of  all  identified  IR-active  species  is  given  in  Table  1, 
and  the  temporal  evolution  of  the  mole  abundances  of  gaseous 
species  is  shown  in  Fig.  3.  N02,  H20,  N20,  and  NO  are  quantified 
based  on  the  H1TRAN  database  [31];  HONO,  CH3N3,  CH30N02  and 
CH3ONO  are  quantified  based  on  data  in  the  literature  [32-35] 
on  the  gas-phase  IR  spectra  of  these  species;  CH3N02  is  quantified 
in  this  work  by  obtaining  a  reference  IR  spectrum  after  evaporating 
0.2  pi  liquid  CH3N02  (purchased  from  Sigma-AIdrich)  into  the 
same  chamber  which  is  introduced  in  Fig.  1,  and  the  mole  fraction 
of  N204  is  estimated  based  on  the  assumption  that  equilibrium  is 
established  quickly  between  N02  and  N204  [26].  The  aerosol  prod¬ 
uct,  monomethylhydrazinium  nitrite,  has  not  been  quantified  due 
to  lack  of  reference  calibration  data.  It  should  be  noted  that  any 
significant  amount  of  MMH  has  not  been  detected  in  any  of  the 
390  spectra  acquired,  which  suggests  that  either  the  chemical  reac¬ 
tion  rate  of  MMH  with  N02/N204  is  much  faster  than  the  evapora¬ 
tion  rate  of  MMH  liquid,  or  that  the  tip  of  the  needle  on  the  syringe 
used  for  injection  does  not  penetrate  any  portion  of  the  volume  of 
the  modulated  beam  of  the  FTIR  spectrometer,  or  both.  By  carefully 
examining  the  acquired  spectra  shown  in  Fig.  2b-d,  the  Q-branch 
of  several  bands  attributed  to  MMH  can  be  detected,  however. 

Figure  4a-c  shows  selected  FT-1R  spectra  obtained  from  a  typi¬ 
cal  test  with  the  chamber  initially  filled  with  5.8%  (molar  basis) 
MMH  vapor  and  94.2%  N2.  The  spectrum  of  MMH  acquired  before 
introducing  N02/N204  is  shown  in  Fig.  4a  and  detailed  assignments 
of  MMH  vibrations  are  provided  in  previous  work  [15].  Figure  4b 
shows  the  spectrum  acquired  at  the  moment  when  N02/N204/N2 
mixture  is  injected  into  the  chamber,  and  Fig.  4c  shows  the  last 
spectrum  acquired  in  the  test.  Figure  4d  is  acquired  by  subtracting 
the  IR  spectral  bands  of  MMH  and  N02  from  Fig.  4c.  The  dominant 
species  identified  in  Fig.  4d  include  monomethylhydrazinium  ni¬ 
trite  and  methyl  diazene  (CH3N=NH).  Very  small  amounts  of 
N20,  CH3N3  and  CH3ON02  were  also  identified.  A  summary  of  these 
identified  species  and  their  major  IR  frequencies  are  given  in 
Table  1. 

It  is  generally  believed  that  the  first  step  in  reactions  between 
MMH  and  N02  is  direct  H-atom  abstraction  by  N02  leading  to  for¬ 
mation  of  HONO  and  a  methylhydrazyl  radical  (CH3NHNH  or 
CH3NNH2)  [13,36],  This  is  supported  by  our  observation  of  HONO 
formation: 

CH3NHNH2  +  N02  CH3NNH2  (or  CH3NHNH)  +  HONO  (1 ) 

The  methylhydrazyl  radical  can  further  react  with  N02  with  poten¬ 
tial  pathways  discussed  by  Catoire  et  al.  [13]  and  Ishikawa  and 
McQuaid  [9],  Catoire  et  al.  suggested  that  CH3N(N02)NH2  and 
CH3N(ONO)NH2  are  the  major  stable  pre-ignition  products  which 
account  for  the  formation  of  a  white-yellow  condensate  [43].  How¬ 
ever,  these  species  have  not  been  identified  as  major  early  species  in 
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Fig.  2.  Selected  IR  spectra  from  a  typical  oxidizer-rich  test  with  N02/N204:  (a)  IR  spectrum  of  N02/N204  mixture  acquired  before  injecting  MMH,  (b)  IR  spectrum  acquired  at 
t  =  7.5  s  after  injecting  MMH,  (c)  IR  spectrum  acquired  at  t  =  15  s  after  injecting  MMH,  (d)  IR  spectrum  obtained  by  subtracting  H20,  N02,  N204,  N20  and  CH30N02  from  the  IR 
spectrum  in  (c),  and  (e)  IR  spectrum  of  monomethylhydrazinium  nitrate  acquired  in  an  early  work  [15]. 


Table  1 

IR  frequencies  of  species  discussed  in  this  work  (vs  =  very  strong,  s  =  strong,  m  =  medium). 


Species 

Major  IR  bands  observed  in  this  work  (cm  *) 

Refs. 

Reactants: 

MMH 

3261  (s),  2960(vs),  2852(vs),  2785(vs),  1595(m),  1480(s),  1451(s),  1293(m),  1119(s),  968(s),  888(vs),  771(vs) 

[37] 

no2 

2907(m),  1617(vs) 

[38] 

n2o4 

1750(s),  1261(s),  749(s) 

[39] 

Products: 

MMH  Nitrite 

3310(s),  3180(s),  3050(s),  2778(s),  2495(m),  1325(s),  1213(vs),  990(m) 

trans-HONO 

3590(m),  1699(s),  1264(s),  791(s) 

[40] 

cis-HONO 

1640(s),  853(s) 

[40] 

CH3N=NH 

3125(m),  2988(m),  2924(m),  1457(m),  1140(m),  844(s) 

[41] 

ch3n3 

2190(m),  2107(s),  1276(m) 

[33] 

ch3ono2 

2965(m),  1666(vs),  1435(m),  1290(vs),  1017(s),  855(s),  759(m),  660(m) 

[34] 

ch3ono 

1678(s),  1620(vs),  991(s),  811(vs) 

[35] 

ch3no2 

1 583(s),  1397(s),  1378(s) 

[42] 

h2o,  no,  n2o 

See  HITRAN  database 

[38] 
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Fig.  3.  Temporal  evolution  of  the  mole  abundances  of  gaseous  species  when  0.2  pi  MMH  was  injected  into  the  chamber  reactor  and  reacted  with  N02/N204  in  1  atm  N2  (Note: 
the  aerosol  product,  MMH-HONO,  is  not  shown  in  this  plot  since  it  has  not  been  quantified  in  this  work). 


Fig.  4.  Selected  1R  spectra  from  a  typical  test  with  MMH  injected  first;  that  is,  the  fuel-rich  case:  (a)  IR  spectrum  of  MMH  vapor  acquired  before  introducing  the  oxidizer,  (b)  1R 
spectrum  acquired  at  the  moment  of  introducing  NO2/N2O4/N2  mixture,  (c)  IR  spectrum  acquired  at  t  =  1 9.5  s  and  (d)  1R  spectrum  obtained  by  subtracting  MMH  and  N02  from 
the  IR  spectrum  in  Fig.  4c. 


the  present  work.  Based  on  thermochemical  calculations,  Ishikawa 
and  McQuaid  [9]  suggested  alternative  pathways  of  direct  H  atom 


abstraction  from  the  methylhydrazyl  radical  by  N02  leading  to 
the  formation  of  methyl  diazene  (CH3N=NH)  or  CH2=NNH2.  We  ob- 
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serve  CH3N=NH  only  as  a  major  product  (Fig.  4d),  and  therefore  we 
conclude  that  the  likely  reaction  between  methylhydrazyl  radical 
and  N02  is  HONO  elimination  with  the  formation  of  methyl  diazene 
(CH3N=NH). 

CH3NNH2  (or  CH3NHNH)  +  N02  —  CH3N  =  NH  +  HONO  (2) 

It  should  be  noted  that  with  excess  N02,  CH3N=NH  can  be  rap¬ 
idly  oxidized  [13],  and  is  not  detected  in  Fig.  2. 

3.2.  Aerosol  formation 

In  previous  studies  [16,43],  a  large  amount  of  condensate  was 
observed  when  mixing  gaseous  MMH  and  N02.  The  condensate 
was  collected  and  examined  ex  situ  with  IR  spectroscopic  tech¬ 
niques  to  reveal  it  contained  monomethylhydrazinium  nitrate  in 
addition  to  other  species.  Similar  findings  were  observed  in  reac¬ 
tions  between  liquid-phase  MMH  and  gaseous  N02  by  Mayer 
et  al.  [17]  In  liquid-phase  reactions  between  MMH  and  N204  di¬ 
luted  by  CCI4,  Saad  et  al.  identified  many  different  products,  includ¬ 
ing  methanol,  methylamine,  dimethylamine,  methylnitrosamine, 
dimethylnitrosamine  and  N-methylformamide  in  addition  to 
MMH  HN03  [18].  These  findings  from  Saad  et  al.  suggest  that  the 
reaction  pathways  in  the  liquid  phase  are  quite  different  than  those 
in  the  gas  phase,  as  indicated  by  the  species  identification  in  Figs.  2 
and  4.  Similar  to  the  results  by  Seamans  et  al.  [16],  we  also  observe 
an  in  situ  formation  of  a  condensate  (fog)  when  MMH  vapor  reacts 
with  N02/N204/N2  mixture.  This  observation  is  from  a  separate  test 
in  which  the  two  reactants  were  introduced  into  a  transparent 
plastic  chamber.  An  in  situ  condensate  formation  can  also  be  in¬ 
ferred  from  changes  in  the  baseline  slope  in  the  IR  spectra  in 
Fig.  2,  and  Kravets  et  al.  [44]  present  an  excellent  discussion  on 
the  effect  of  small  particles  on  spectral  transmittance.  The  IR  trans¬ 
mittance  decreases  with  increasing  wavenumber  due  to  the  Ray¬ 
leigh  scattering  of  small  particles  which  generates  a  baseline 
slope  in  the  IR  spectra.A  preliminary  analysis  of  the  condensate 
by  Seamans  et  al.  [16]  indicated  a  mass  composition  of  C:  14.6%, 
H:  8.1%,  N:  46.4%,  and  0:  30.9%  (by  difference).  Catoire  et  al.  [13] 
suggested  that  the  monomethylhydrazinium  nitrate  is  not  likely 
formed  in  gas-phase  reactions.  Instead  Catoire  et  al.  attributed  this 
condensate  to  covalently  bonded  nitro  and  nitrite  compounds  pro¬ 
duced  from  methylhydrazyl  radical  and  N02  because  these  species 
have  calculated  a  mass  composition  of  C:  13.2;  H:  5.5;  N:  46.1  and 
0:  35.1,  consistent  with  the  analysis  of  Seamans  et  al.  [16]  Our 
in  situ  IR  analysis  reveals  that  this  condensate  is  mainly  composed 
of  an  ionic  monomethylhydrazinium  nitrite  salt  (MMHHONO). 
This  nitrite  salt  has  a  calculated  mass  composition  of  C:  12.9;  H: 
7.5;  N:  45.2  and  0:  34.4,  also  in  agreement  with  Seaman’s  analysis. 
Based  on  HONO  as  a  weak  acid  [45]  (p/Ca  =  3.16)  and  MMH  as  a 
well-known  base  [46]  (p/Ca  =  8.1),  formation  of  monomethylhydr¬ 
azinium  nitrite  through  an  acid-base  proton  transfer  reaction  is 
very  likely.  The  gas-phase  acid-base  reaction  has  been  thoroughly 
investigated  by  atmospheric  chemists,  [47-50]  and  it  is  generally 
believed  that  a  hydrogen-bonded  neutral  complex  is  formed  prior 
to  the  occurrence  of  proton  transfer.  This  type  of  reaction  is  widely 
known  as  ‘gas-to-aerosoT  or  ‘gas-to-particle’  reaction  because  an 
aerosol  of  salts  is  usually  formed. 

CH3NHNH2  +  HONO  -►  CH3NHNH2  ■  ■  ■  H  ■  ■  ■  ONO 

->  [CH3NHNH3  •  ONO“]  (3) 

[CH3NHNHJ  •  0N0~]n  ,  +  [CH3NHNH3  •  ONO-] 

-*■  [CH3NHNH+  ■  0N0~]n  (4) 

Seamans  et  al.  [16]  found  that  the  amount  of  such  condensate 
was  more  than  20%  by  weight  of  the  total  reactants  in  some  exper¬ 
iments.  Therefore,  reactions  (3)  and  (4)  must  be  included  in  the 


MMH/RFNA  mechanism.  Furthermore,  the  condensate  formation 
can  play  an  important  role  in  the  ignition  event  for  several  reasons. 
First,  the  associated  proton  transfer  is  quite  exothermic,  thus  pro¬ 
viding  energy  for  increase  in  temperature  of  the  system  and  energy 
for  the  endothermic  H-abstraction  reactions  forming  HONO.  Sec¬ 
ond,  we  consider  that  the  formed  particle  provide  a  suitable  med¬ 
ium  for  enhanced  surface  reactions.  Third  and  finally,  as  molecules 
become  “trapped”  in  a  condensed  phase,  the  concentration  of  gas¬ 
eous  species,  such  as  MMH,  N02  and  HONO,  is  likely  to  be  in¬ 
creased  leading  to  increased  reaction  rates. 

3.3.  Comparison  between  fuel-rich  and  oxidizer-rich  conditions 

Under  fuel-rich  conditions,  shown  in  Fig.  4,  reactions  (l)-(4)  are 
relevant.  With  excess  MMH,  HONO  formed  by  reaction  (1)  and  (2) 
can  be  almost  fully  consumed  by  reaction  (3)  to  form  MMH  HONO. 
The  amount  of  gaseous  HONO  in  all  spectra  is  almost  unnoticeable 
due  to  the  high  concentration  of  MMH  driving  the  salt  formation 
with  HONO.  Methyl  diazene  formed  by  reaction  (2)  is  one  of  the  fi¬ 
nal  products  with  excess  MMH  shown  in  Fig.  4. 

Under  oxidizer-rich  conditions,  methyl  diazene  (CH3N=NH)  is 
further  oxidized  and  barely  detected,  as  indicated  in  Fig.  2.  Based 
on  IR  analysis  subsequent  reactions  between  CH3N=NH  and  N02 
produces  large  amounts  of  CH3ONO  and  CH3N02.  The  first  step  is 
the  likely  H-atom  abstraction  by  N02  leading  to  the  formation  of 
methyl  diazene  radical  CH3N=N  and  HONO.  The  possible  reactions 
between  CH3N=NH  and  N02  under  oxidizer-rich  conditions  are: 

CH3N=NH  +  N02  CH3N=N'  +  HONO  (5) 

CH3N=N'  ->  CHj  +  N2  (6) 

CHj  +  N02  -»  CH3NO2  (or  CH3ONO)  (7) 


4.  Theoretical  results 

4.1.  H-atom  abstraction  from  MMH 

The  oxidation  of  MMH  via  sequential  HONO  formation  and  the 
final  N2  generation  were  studied  computationally  as  shown  in 
Scheme  1.  The  first  H-abstraction  from  MMH  has  been  studied  in 
detail  by  McQuaid  and  Ishikawa  [7],  Their  reported  barrier  heights 
for  H-abstractions  from  three  different  positions  (H  on  N-CH3  and 
two  H  on  N-NCH3  -  cis  and  trans  to  the  methyl  group,  respectively) 
to  form  cis-HONO  are  10.1, 10.6  and  11.2  kcal/mol  (without  ZPE  or 
thermocorrections)  at  CCSD(T)/6-311+G(2df,p)//CCSD/6-31+G(d,p) 
level  relative  to  free  N02  and  MMH. 

These  results  are  close  to  outvalues  of  8.0,  9.8  and  10.3  kcal/mol 
(including  ZPE  and  temperature  corrections,  relative  to  free  N02 
and  MMH).  We  find  that  the  most  readily  abstracted  H-atom  is  from 
the  nitrogen  with  the  methyl  group,  in  agreement  with  the  previous 
study.  The  binding  energy  between  the  product  HONO  and  the 
free  radical  is  substantial  (10-13  kcal/mol  enthalpy),  reducing  the 
free  energy  of  the  post-reaction  complex  to  even  lower  levels  than 
the  unbound  free  radicals  Glntl  -  Glnt3.  However  in  normal  exper¬ 
imental  conditions  the  partial  pressure  of  HONO  is  usually  much 
lower  than  1  atm  (the  reference  state),  favoring  the  formation  of 
unbound  free  radicals,  after  taking  concentration  correction  into  ac¬ 
count.  The  interconversion  between  GInt2  and  Gint3  via  N— N  bond 
rotation  has  a  high  barrier  (~29  kcal/mol)  because  the  N— N  bond 
has  some  double-bond  character  due  to  the  delocalization  of  the 
N  lone  pair.  Abstraction  of  an  H  atom  from  a  methyl  group  (GTS4) 
is  ~10  kcal/mol  higher  in  energy  than  from  N  atoms.  Furthermore, 
this  reaction  is  16.1  kcal/mol  endothermic,  in  contrast  to  H-abstrac¬ 
tions  from  N  atoms,  which  are  almost  thermoneutral.  Therefore 
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Scheme  1.  Reactions  between  MMH  and  N02  in  gas  phase.  Barriers  of  all  H-abstractions  from  N  of  MMH  to  form  HONO  are  about  10  keal/mol,  and  the  same  H-abstraction 
from  methyl  group  is  16.1  keal/mol  endothermic  and  has  -  1 0  keal/mol  higher  barrier,  rendering  the  oxidation  of  carbon  slower  at  low  temperature.  Enthalpies  and  Gibbs  free 
energies  (in  the  parentheses)  are  calculated  at  298.15  K  and  1  atm.  cis-HONO  at  standard  state  is  used  as  reference  product. 


H-atom  abstraction  from  the  methyl  group  will  not  play  an 
important  role  during  the  pre-ignition  event,  and  experimentally 
we  observed  abundant  methyl-containing  compounds,  indicating 
the  inertness  of  methyl  group  at  room  temperature.  Consequently 
reactions  beyond  Glnt4  are  not  considered  in  Scheme  1. 


4.2.  H-atom  abstraction  from  CH3NNH2  or  CH3NHNH 

The  H-abstraction  reactions  from  either  CH3NNH2  or  CH3NHNH 
very  likely  have  no  barrier  because  of  stabilization  of  the  N  p  orbi¬ 
tal  bonded  to  the  H  with  delocalization  of  the  adjacent  N  lone  pair. 
Indeed  we  could  not  find  a  transition  state  in  the  electronic  energy 
surface  despite  an  exhaustive  search.  The  free  energy  surface  after 
including  ZPE  might  well  lead  to  a  barrier.  Ishikawa  and  McQuaid 
[9],  using  the  MPWB1  I<  functional,  found  that  1-2  keal/mol  of  ki¬ 
netic  energy  is  enough  to  activate  H-atom  abstraction,  indicating 
very  low  barriers  for  H-abstraction.  GInt7  -  GlntlO  are  products 
from  recombination  between  two  radicals  (Glntl  -  Glnt3  and 
N02).  All  of  them  are  enthalpically  and  entropically  less  stable  than 
the  product  of  HONO  formation,  GInt5  and  GInt6.  These  recombi¬ 
nation  products  are  not  detected  in  our  experiment,  where  MMH 
and  N02  are  in  gas  phase  diluted  with  N2„  With  higher  concentra¬ 
tion  of  MMH  and  N02  or  reacting  in  liquid  phase  will  change  free 
energy  profile,  making  the  recombination  products  more  favorable. 


4.3.  H-atom  abstraction  from  CH3N=NH 

The  enthalpic  barrier  for  the  H-atom  abstraction  is  11.1  kcal/ 
mol  for  anti-CH3N=NH  (Glnt5  ->  GTS6  ->  Glntl  1 )  and  6.2  kcal/ 
mol  for  syn-CH3N=NH  (GInt6  -*  GTS7  -»  Glntl  1).  The  CH3N=N 
radical  (Glntl  1)  can  either  break  the  C— N  bond  to  release  N2  and 
CH3  radical  via  TS9  (ARL  mechanism  reaction  No.  456)  with  only 
2.4  keal/mol  barrier,  or  undergo  N02  attack  at  different  orienta¬ 
tions  to  form  CH3N02  and  CH3ONO  via  TS8  and  TS10  with  barrier 
heights  2.5  and  3.1  keal/mol,  respectively. 

It  is  not  likely  that  the  free  radicals  CH3  and  CH3NN  play  a 
significant  role  in  the  formation  of  CH30N02.  N02  is  much  more 
abundant  than  ONON02  in  the  gas  phase  (see  Section  4.5),  there¬ 
fore  the  recombination  of  CH3  or  CH3NN  with  N02  to  give  CH3N02 
or  CH3ONO  should  be  much  faster  than  the  recombination  of  CH3 
or  CH3NN  with  ONON02.  The  other  pathway  to  form  CH30N02  is 
via  the  recombination  of  N02  with  the  CH30  radical,  which  is  gen¬ 
erated  from  the  reaction  CH3+N02  ->  CH30+NO;  however,  this 
reaction  cannot  fully  explain  the  production  of  CH3ON02  either 
due  to  the  much  different  concentrations  of  NO  and  CH3ON02.  In 
addition,  the  formation  of  CH2N2  via  H-abstraction  from  methyl 
group  is  also  found  to  have  a  low  barrier  (3.3  keal/mol).  However 
this  product  was  not  detected  experimentally,  indicating  that  the 
direct  dissociation  of  methyl  free  radical  from  Glntl  1  may  be  much 
faster  than  other  bimolecular  pathways. 


W.-G.  Liu  et  al. /  Combustion  and  Flame  1 60  (2013)  970-981 


977 


4.4.  Formation  of  MMHHONO  aerosol  and  its  lower  reactivity 

In  each  step  of  Scheme  1,  HONO  is  produced,  which  is  able  to 
form  aerosol  with  unreacted  MMH  because  the  basic  N  atoms  on 
MMH  are  able  to  accept  a  proton  from  HONO.  As  shown  in 
Scheme  2,  protons  in  MMH  HONO  complexes  in  gas  phase  prefer 
to  stay  on  HONO.  The  enthalpies  to  form  nHONO  MMH  complexes 
(n  =  1, 2)  are  roughly  additive  (-21.5  kcal/mol  from  the  sum  of  en¬ 
thalpy  of  GIntl2  and  13  vs.  -18.9  kcal/mol  of  GIntl4)  and  the  free 
energies  to  form  these  complexes  are  about  thermoneutral  (-2.1, 
0.1  and  0.7  kcal/mol  for  Glntl2,  GIntl3  and  GIntl4).  In  the  strong 
solvation  environment  (in  the  aerosol),  the  proton  transfer  from 
HONO  to  MMH  is  more  favored  and  exothermic,  which  accounts 
for  the  observation  of  ONO~  anion  in  the  IR  spectra.  However, 
the  reaction  heat  to  transfer  one  and  two  protons  from  the 
MMH-2HONO  complex  in  a  solvated  system  differ  little  (5.7  vs. 
8.0  kcal/mol  exothermic),  unlike  complexes  in  the  gas  phase.  It  is 
because  in  solvated  system  the  electrostatic  repulsion  between 
two  positively  charged  N  centers  partially  cancels  the  energy  gain 
from  the  neutralization.  Hence,  the  pathways  for  continued  growth 
and  composition  of  the  aerosol  can  be  either  MMH2+-20NCT  or 
MMH+ONO~.  In  the  oxidizer-rich  case,  however,  the  former  com¬ 
position  is  more  likely.  The  spontaneous  nucleation  in  gas  phase 
followed  by  the  exothermic  growth  make  the  aerosol  formation  a 
rapid  process  as  observed  in  the  experiment. 

We  also  found  that  the  MMH  HONO  aerosol  is  less  reactive  than 
free  MMH.  Reactions  to  abstract  H  from  MMH  HONO  complex  have 
barriers  7-9  kcal/mol  higher  than  the  same  H-abstraction  from 


MMH.  The  reason  is  that  in  free  MMH,  the  lone  pairs  on  the  N 
can  stabilize  the  transition  state  via  resonance,  whereas  such  sta¬ 
bilization  is  less  available  when  lone  pairs  donate  electron  density 
to  the  proton,  resulting  in  a  higher  barrier.  The  increase  of  the  bar¬ 
rier  for  HONO  formation  resulting  from  the  salt  formation  was  also 
observed  in  the  case  of  alkylamine  and  nitric  acid  [51].  As  a  result, 
growth  of  the  particles  is  favored  over  H-abstraction  reactions  at 
the  low  temperatures  of  our  experiments.  As  the  temperature  in¬ 
creases,  particles  will  either  undergo  H-abstraction  reactions  or 
evaporate  into  MMH  and  HONO  so  that  particles  will  shrink  and 
eventually  disappear.  Such  a  phenomenon  of  particle  disappear¬ 
ance  is  observed  just  prior  to  ignition  in  the  MMH  HN03  reaction 
system  [15], 

4.5.  Reactions  of  asymmetric  dimer  of  N02,  0N0N02  in  gas  phase 

The  formation  of  0N02  and  CH30N02  cannot  be  explained  by 
simple  H-atom  abstractions  or  recombination  of  the  radical  inter¬ 
mediates  with  N02.  One  potential  source  of  nitrate  is  from  the 
isomerization  of  N204.  It  is  known  that  liquid  N02  dimerizes  to 
form  N204  and  disproportionates  into  NO+  and  NOJ  [52],  N02 
can  also  react  with  water  vapor  to  give  HONO  and  HN03.  Finlay- 
son-Pitts  and  coworkers  [53]  proposed  that  the  asymmetric  isomer 
0N0N02  is  the  key  intermediate  as  the  source  of  nitrate.  Our  pre¬ 
vious  study  shows  that  the  reaction  to  form  0N0N02  has  low 
enthalpic  barrier  (<5  kcal/mol)  [54],  which  means  equilibrium  be¬ 
tween  N02  and  0N0N02  is  very  fast.  Recently  Lai  et  al.  [55]  also 
found  that  0N0N02  can  play  an  important  role  in  the  hypergolic 
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reaction  between  hydrazine  and  liquid  NTO.  We  would  also  like  to 
check  if  ONON02  plays  the  similar  role  in  the  gas  phase  reaction 
between  MMH  and  N02. 

The  easiest  reactions  between  MMH  and  trans-ONON02  is  the 
new  N— N  bond  formation  between  NO+  and  electron-rich  N  atoms 
on  MMH  followed  by  the  proton  transfer  from  N— H  bond  to  N03  to 
form  nitric  acid,  as  shown  in  Scheme  3.  Although  it  is  very  easy  to 
form  nitric  acid  (enthalpic  barrier  is  about  1  kcal/mol),  to  form 


methyl  nitrate  is  difficult  in  gas  phase  (34  kcal/mol  enthalpic  bar¬ 
rier)  due  to  the  unfavorable  charge  separation  (nitrate  anion  has  to 
attack  from  the  back  of  methyl  group,  which  is  far  from  the  posi¬ 
tive  N  center).  Therefore,  the  gas  phase  reactions  between  MMH 
and  ONON02  can  easily  generate  HN03  and  explain  the  experimen¬ 
tal  observation  of  N03 ,  however,  they  cannot  explain  the  abundant 
CH30N02  observed  in  the  IR  spectrum,  which  we  surmised  to  be 
produced  on  the  surface  of  aerosol,  as  studied  in  the  Section  4.6. 
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Scheme  4.  The  reactions  between  ONON02  and  MMH  in  water  to  simulate  potential  energy  surface  in  the  aerosol.  Enthalpies  are  calculated  at  298.15  K  and  1  atm.  Several 
gas  products,  such  as  CH30N02  and  N20,  can  desorb  from  the  aerosol  and  be  observed  via  IR  spectra. 
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Scheme  6.  The  reactions  between  CH3N3  and  0N0N02  in  gas  phase.  All  barriers  are  significantly  higher  than  the  corresponding  reactions  in  water,  indicating  the  importance 
of  solvation  effect.  Enthalpies  and  Gibbs  free  energies  (in  the  parentheses)  are  calculated  at  298.15 1<  and  1  atm. 


0N0N02  can  also  react  with  intermediates  from  HONO  formation, 
such  as  GInt5  and  GInt6,  as  shown  in  Scheme  3.  N  atoms  on  GInt5 
and  Glnt6  are  sp 2  hybridized,  not  as  electron-rich  as  the  sp3  hybrid¬ 
ized  ones,  therefore  barriers  to  abstract  H  from  Glnt5  and  Glnt6  to 
form  HN03  are  about  10  kcal/mol  higher  than  H-abstractions  from 
MMH. 

4.6.  Reactions  facilitated  by  aerosol  to  form  CH30N02 

Several  experiments  indicate  that  N02  can  react  with  water  or 
alcohol  heterogeneously  (surface-catalyzed)  [53,56-58],  probably 
via  the  same  0N0N02  intermediate  [53].  The  aerosol  of  ONCT 
and  MMH  cation  (such  as  Sintl-Sint4)  provides  polar  ionic  surface 


that  can  stabilize  the  asymmetric  0N0N02  and  promote  nitrate 
formation.  This  pathway  can  be  favored  at  higher  N02  concentra¬ 
tions  because  0N0N02  concentration  is  proportional  to  [N02]2, 
which  also  explain  the  experimental  observation  that  the  forma¬ 
tion  of  methyl  nitrate  was  only  observed  in  the  N02-rich  atom- 
sphere.  Because  asymmetric  0N0N02  has  high  dipole  moments 
(3.45D  for  cis  and  2.96D  for  trans  at  B3LYP/6311G*+  level)  [54] 
we  expect  the  concentration  of  0N0N02  to  be  greatly  increased 
if  it  is  absorbed  on  a  polar  surface  or  formed  a  molecular  complex 
with  a  polar  molecule. 

To  study  reactions  between  0N0N02  and  MMH  on  the  polar 
surface  or  the  sub-surface  of  the  aerosol,  we  used  an  implicit  solva¬ 
tion  model  to  include  the  interaction  between  reactants  and 
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surrounding  ions.  Our  proposed  reaction  mechanism  between 
MMH  and  0N0N02  and  the  corresponding  enthalpies  are  shown 
in  Scheme  4. 

In  solution  phase,  0N0N02  has  strong  tendency  to  dissociate 
into  NO3  and  NO+  with  the  presence  of  electron-rich  MMH.  The 
electron-deficient  NO+  makes  N— N  bond  with  the  electron-rich 
N,  preferably  the  methyl  substituted  N,  on  MMH  to  form  SInt5 
and  Slnt6.  From  Slnt5,  nitrate  anion  can  attack  from  the  back 
methyl  group  to  form  methyl  nitrate  with  35.7  kcal/mol  barrier. 

Once  the  N— N  bond  is  formed,  the  acidity  of  the  adjacent  N-H 
bond  increases  and  the  overall  complex  is  similar  to  the  salt  gener¬ 
ated  from  the  neutralization  between  HN03  and  corresponding 
amine.  We  assumed  that  the  proton  transfer  in  such  a  polar  envi¬ 
ronment  has  a  low  barrier.  Therefore  to  study  tautomers  as  inter¬ 
mediates  is  enough  to  depict  the  potential  energy  surface.  In  our 
calculations,  these  tautomers  are  tightly  constrained  to  conserve 
the  number  of  protons,  while  in  aerosol  protons  can  exchange  with 
the  environment,  implying  our  barriers  should  be  the  upper 
bounds  of  true  barriers. 

Even  with  free  nitrate  anion  in  aerosol,  the  barriers  for  the 
nucleophilic  substitution  of  the  methyl  group  on  MMH  are  consid¬ 
erable.  Among  several  possible  tautomers,  Slnt8  and  Slntl  1  are  the 
easiest  two  to  form  CH30N02  via  STS3  and  STS4  to  SInt9,  with  bar¬ 
riers  36.2  and  33.3  kcal/mol,  respectively,  giving  the  final  products, 
hydrazoic  acid  and  water.  See  SI  for  other  tautomers  and  corre¬ 
sponding  TS  to  form  CH30N02.  It  is  also  possible  to  form  N20 
and  methyl  amine  by  breaking  the  N— N  bond  in  Slntl  2  via  STS5 
with  13.7  kcal/mol  barrier,  corresponding  to  the  experimental 
observation  of  N20. 

An  easier  path  to  form  CH30N02  involves  the  intermediate 
CH3N3  formed  via  STS6  with  15.9  kcal/mol  enthalpic  barrier.  As 
shown  in  Scheme  5,  in  solution  phase,  CH3N3  prefers  to  react  with 
0N0N02  to  form  a  5-member  ring  intermediate  (SIntl5)  with  neg¬ 
ligible  barrier,  which  then  decomposes  to  release  N2  with 
14.8  kcal/mol  barrier,  and  the  product  (Slntl 6)  undergoes  nucle- 
ophlic  attack  of  nitrate  to  form  CH30N02  and  N20  with  17.0  kcal / 
mol  enthalpic  barrier  and  32.3  kcal/mol  exothermicity.  This  mech¬ 
anism  fits  to  the  fact  that  azide  reacts  completely  with  NOC1  or 
N204  in  condense  phase  to  give  nitrate,  N20  and  N2  [59,60],  and 
was  proposed  in  Ref.  [18].  The  5-member  ring  intermediate  and 
its  decomposition  pathway  have  also  been  reported  [61], 

The  solvation  effect  plays  an  important  role  in  facilitating 
the  decomposition  of  CH3N3,  because  in  gas  phase,  similar  reac¬ 
tions  to  produce  CH30N02  have  enthalpic  barriers  higher  than 
25  kcal/mol,  as  shown  in  scheme  6.  In  gas  phase,  to  form  the  same 
5-member  ring  intermediate  via  GTS22  has  a  27.5  kcal/mol  enthal¬ 
pic  barrier,  The  STS1 1-like  transition  state,  GTS24,  does  not  con¬ 
nect  to  a  5-member  ring  intermediate  in  the  MEP  scan,  instead  it 
leads  to  the  path  to  dissociate  N2  from  CH3N3  directly,  with  barrier 
height  29.1  kcal/mol.  These  reaction  paths  do  not  lead  to  rapid  pro¬ 
duction  of  CH30N02,  leading  to  the  conclusion  that  the  formation 
of  CH3N3  and  CH30N02  take  place  on  the  surface  of  aerosol  but 
not  in  gas  phase. 

Once  N20,  CH30N02,  and  CH3N3  are  formed  on  the  surface  of 
aerosol,  they  can  desorb  and  drive  reactions  further  towards  com¬ 
pletion.  The  energetics  is  favorable  due  to  the  partial  oxidation  of 
MMH.  Also  the  exchange  between  gas  phase  and  aerosol  species 
is  plausible  -  N02  can  be  absorbed  onto  the  aerosol  surface  and 
oxidize  intermediates  shown  in  Scheme  2,  leading  to  a  complicated 
multiphase  picture  of  pre-ignition  reactions. 


5.  Conclusions 

We  studied  reactions  between  MMH  and  N02  vapor  in  a 
gold-coated  chamber  reactor  with  Fourier  transform  infrared 


spectrometry  at  both  MMH-rich  and  N02-rich  conditions.  At  low 
concentration  of  N02  the  major  products  are  MMH  HONO  and 
CH3N=NH  and  the  minor  products  are  N20,  CH3N3,  CH3N02  and 
CH3ONO.  Our  QM  calculations  elucidate  possible  mechanisms  of 
H-abstraction  by  N02  from  MMH  to  form  HONO,  which  then  forms 
condensate  by  reacting  with  MMH  in  fuel-rich  condition.  We  find 
that  CH3N=NH  is  formed  after  a  second  H-abstraction  from 
MMH.  Further  H-abstraction  produces  CH3ONO,  CH3N02  and  N2. 
At  higher  concentrations  of  N02  the  major  products  are  mono- 
methylhydrazinium  nitrite  and  methyl  nitrate.  The  formation  of 
methyl  nitrate  is  attributed  to  the  asymmetric  isomer  0N0N02  of 
N204,  which  is  favored  at  high  N02  concentration.  Our  ab  initio  cal¬ 
culations  indicate  that  further  reactions  between  MMH  and  ONON- 
02  facilitated  by  the  surface  of  the  aerosol  or  inside  can  generate 
CH30N02,  CH3N3  and  N20,  products  observed  in  N02-rich  experi¬ 
mental  conditions.  This  study  illustrates  the  heterogeneous  nature 
of  the  pre-ignition  reactions  between  MMH  and  N02. 
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